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Abstract 

We compute the renormalon ambiguity of the static potential, in the hmit 
of a large number of flavors. An extrapolation of the QED result to QCD 
implies that the large distance behavior of the quark potential is arbitrary 
in perturbation theory, as there are an infinite number of prescriptions to 
assign. The shape of the potential at large distances is not only affected by the 
renormalon pole closest to the origin of the Borel plane, but a resummation 
of all renormalon contributions is required. In particular, confinement can 
be accommodated, but it is not explained. At short distances there is no 
indication of a linear term in the potential. 
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A way to learn about intrinsic limitations of perturbation theory in the determination of 
a physical quantity is via the study of the large order behavior of the perturbative expansion. 
Although the coupling constant depends only logarithmically on the scale, the perturbative 
expansion can produce power suppressed corrections, as it is an asymptotic series. Some 
nonperturbative features of the theory can be inferred from the positions and residues of the 
singularities in the Borel transform, with respect to the coupling constant 

At present it is not known how to compute renormalons in full QCD. Most calculations 
have been performed by summing fermion bubble insertions in the gluon propagator. The 
neglected graphs are presumably just as important as those taken into account, so the 
result of summing the fermion bubble chain does not prove the existence of renormalons in 
QCD. Calculating renormalons via bubble summation is rigorous in the limit of QED with 
a large number of fermions, Nf — oo. In this case a fermion bubble insertion in a photon 
line yields a factor of a = aNf, which can be summed to all orders. The contribution of 
other diagrams is suppressed by powers of l/A^/- Extrapolating the large Nf QED result 
to QCD is accomplished by assuming that the ultraviolet and infrared properties can be 
interchanged. Such studies have received considerable attention in the literature recently, 
mostly in connection with the heavy quark expansion (see, for example IQ-^)- 

In this letter we study renormalon ambiguities in a static property of QCD, i.e., the 
potential between a heavy quark and an anti-quark. Neglecting color factors, the leading 
order one gluon exchange diagram contributes to the potential by the usual Coulomb term 



Diagrams with multiple gluon exchange are suppressed by powers of l/Nf, while the insertion 
of fermion bubbles in the gluon line involves the coupling a, which is treated as a parameter 
of order unity. Therefore, we have to consider an arbitrary number of bubbles inserted into 
the one gluon exchange diagram. 

If the perturbative expansion of the potential is given by 




(1) 
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then the Borel transform is defined as 

nr) = t^u-. (3) 

n=0 

Here u is the Borel parameter. For convenience, we defined the Borel transformation with 
respect to 60 a, where 60 is the first coefficient of the /5-function, (3{a) = fi^da/dfi'^ = 
-bo + 0{a^). In QED 60 = -iV//37r, while in QCD bo = (11 - |iV/)/47r 0. The reason 
for studying the Borel transform V{r) is that it may converge even if the perturbative series 
for V{r) in Eq. (0) is divergent. In that case the inverse Borel transformation 

POO ^ 

V{r)= / dMe-"/(''»"V(r), (4) 



provides a good definition of V{r). However, if Vn grows at least factorially for large n, then 
V{r) can have singularities on the integration contour. In our case these singularities will 
be isolated poles on the positive real axis that arise from the infrared region of Feynman 
diagrams. These so-called infrared renormalons yield ambiguities in the reconstruction of 
V{r) from V{r), as the result depends on the regularization of the pole contributions. 

The Borel transform of the static potential can be calculated directly, using the Borel 
transform of the resummed gluon propagator |Tl|. This yields 

V{r) = 



bo KfiV J (27r)3 (A;2)i+« 



60 r ' T{2u + 1) 

where /i is the renormalization scale, and C is a regularization scheme dependent constant. 
We work in the momentum subtraction scheme in which C = 0. The singularities of V{r) 
are simple poles along the positive real axis at half-integer values of u 

Uk = k + ^, fc = 0,l,2,.... (6) 

The residues of these poles are 



The physical potential is reconstructed by the inverse Borel transform in Eq. 



Vir) = — f / (A r)^" ' '\ , (8) 

^ ' r boJo ^ ' T{2u + 1) ^ ^ 

where A is the pole of the coupling (the dynamically generated scale), and we have used the 
relation AV/i^ = e-V(feoa[/.2]), 

carry out this integration, the contour has to be deformed 
away from the renormalon poles.|^ Using the QCD /3-function coefficient, we obtain the 
ambiguity in the quark potential 

^''<^'^;^£-<^^>^"'(^- 

Here Ck are arbitrary complex numbers of order unity, related to the ambiguity of regularizing 
the l/{u — Uk) poles of V{r) by their principal value plus Ck times the delta function 6{u—Uk)- 

Inserting the QCD /3-function into Eq. (||) is equivalent to replacing the coupling constant 
evaluated at the fixed subtraction scale in the one gluon exchange diagram with the running 
coupling at the scale of the gluon momentum, a;(A;^). Diagrammatically this corresponds to 
dressing the one gluon exchange diagram with self-energy and vertex corrections computed 
in the leading logarithmic approximation. 

The crucial point is that the static potential is an unambiguous physical quantity. It can 
be calculated on the lattice as the expectation value of a Wilson loop. In view of the above 
ambiguity originating from the renormalon calculation, the static potential can only be free 
of ambiguities if nonperturbative effects cancel Eq. (P). It is in this sense how renormalons 
'trace' nonperturbative effects.^] 

There are a number of points to be made regarding our main result in Eq. (|^): 



*These renormalon poles are not present if we include an infrared regulator, like a gluon mass or 
a vacuum expectation value for the (G^^ G^'^) operator. However, we are interested in what can 
be learned strictly from perturbation theory. 

''^There may also be instanton contributions, and nonperturbative effects that do not correspond 
to any local operator. 



a. In previous calculations, powers of some large mass scale or momentum flow suppressed 
the renormalon contributions of the poles further away from the origin of the Borel plane. In 
our case this scale is replaced by 1/r. At large distances A/(l/r) enhances the contributions 
of the poles far from the origin. Therefore, the ambiguity in the static potential at large 
distances is affected by all renormalon poles. 

b. By varying the values of Ck, the renormalon ambiguity at large distances becomes 
arbitrary. Taking for example Ck — 1 yields 

, ^ , sin(Ar) , 
AV{r) , (10) 

while Ck — (—1)*^ yields 

. , sinh(Ar) , , 

AV{r)r A L_ (11) 

In the first case, the renormalon calculation implies that nonperturbative correction to the 
potential is not confining, while it is confining in the latter case. With a proper choice of 
the Cfc's, an arbitrary potential shape can be produced at large distances. For example, the 
Cjt's can reproduce a nonperturbative correction to the potential of the form A ^Jl + (Ar)^, 
i.e, a linear behavior in r for r ^ A~^. 

c. The renormalon ambiguities do not affect the Coulomb part of the potential, only terms 
that are proportional to r"^^. Near r — the renormalon calculation gives no indication of 
a linearly rising correction to the potential. It seems conceivable to us that this is a real 
feature of QCD, namely 'flux tubes' develop only for r > A~^ Our result does not constrain 
the large r behavior: as was shown above, V{r) can even grow linearly for large r. 

d. There is an analogy with non-renormalizable theories, where knowledge of an infinite 
number of counterterms is needed to derive scattering amplitudes from the Lagrangian. In 
the present case, the values of would specify the gluon propagator at large distances. 

These observations imply that the renormalon calculation supplemented with bubble 
summation has a very limited capability in describing long distance QCD. In this sense our 
results question the suitability of resumming fermion bubble chains for phenomenological 
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estimates of nonperturbative effects in real QCD. One could argue that some values of Ck 
imply confinement, but based on the information available from the renormalon calculation, 
such a choice is completely arbitrary, and is not any better motivated than values of Ck that 
suggest that nonperturbative contributions to the static potential need not be confining. In 
this framework, therefore, the Landau pole gives no indication of what type of nonperturba- 
tive effects occur. Beyond the perturbatively calculable result there are an infinite number 
of parameters that need to be specified to describe the long distance dynamics. 



ACKNOWLEDGMENTS 

We thank Guido Martinelli, John Preskill, and especially Mark Wise for discussions. 



6 



REFERENCES 



[1] G. 't Hooft, in: The Whys of Subnuclear Physics, Proceedings of the 15th International 
School on Subnuclear Physics, Erice, Sicily, 1977, Edited by A. Zichichi (Plenum Press, 
New York, 1979). 

G. Parisi, in: Hadron structure and lepton-hadron interactions, Cargese 1977, Edited 
by M. Levy et al. (Plenum Press, New York, 1979); 
[2] I.l. Bigi, M.A. Shifman, N.G. Uraltsev, and A.I. Vainshtein, TP1-MINN-94/4-T, |hep^ 
ph/9402360| . 



[3] M. Beneke and V.M. Braun, Nucl. Phys. B426 (1994) 301. 
[4] A.V. Manohar and M.B.Wise, Phys. Lett. B344 (1995) 407. 
[5] M.Beneke, V.M. Braun, and V.I. Zakharov, Phys. Rev. Lett. 73 (1994) 3058. 
M. Luke, A.V. Manohar, and M.J. Savage, UTPT 94-21, |hep-ph/ 940 7407 . 



[7] M. Neubert and C.T. Sachrajda, CERN-TH-7312-94, |hep-ph/ 940 7394 
[8] U. Aglietti, GALT-68-1981, |hep-ph/95023'gq . 

[9] L. Maiani, G. Martinelh, and C.T. Sachrajda, Nucl. Phys. B368 (1992) 281; G. Martinelh 
and C.T. Sachrajda, CERN-TH.7517/94, ^ep-ph/ 950253^ . 



[10] H.D. Politzer, Phys. Rev. Lett. 30 (1973) 1346; D.J. Gross and F. Wilczek, Phys. Rev. 

Lett. 30 (1973) 1343. 
[11] M. Beneke, Nucl. Phys. B405 (1993) 424. 



7 



